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1.0  INTRODUCTION 


The  purpose  of  this  study  is  a preliminary  assessment  of  the 
utility  of  the  Nimbus  6 Electrically  Scanning  Microwave  Radiometer 
(ESMR)  for  sea-surface  wind  speed  determination.  The  37  GHz  Nimbus  6 
ESMR,  successor  to  the  19.35  GHz  Nimbus  5 ESMR,  differs  from  the  latter 
not  only  in  frequency  but  also  in  scanning  geometry  and  polarization. 
Previous  studies  [1,  2]  proved  the  utility  of  the  Nimbus  5 ESMR  to  sea- 
surface  wind  determinati jn,  but  also  pointed  out  the  deleterious  effects 
of  a cloudy  atmosphere  and  of  the  not-well -defined  relationship  between 
wind  and  wave  spectrum  on  the  ultimate  estimation  of  wind  from  bright- 
ness temperatures.  The  Nimbus  6 ESMR  with  its  dual  polarization  would 
partially  compensate  for  atmospheric  effects  by  providing  additional 
information  in  the  vertically  polarized  channel.  Nevertheless,  in 
precipitation  areas  and  in  areas  o t liquid  clouds  the  uncertainties 
introduced  by  the  atmosphere  are  still  too  large,  thus  limiting  the 
ability  of  the  Nimbus  6 ESMR  for  accurate  sea-surface  wind  estimation  to 
cloud-free  areas. 

The  utility  of  the  Nimbus  6 ESMR  for  sea-surface  wind  determination 
is  explored  in  this  study  by: 


1.  Analysis  of  wind  equations  derived  from  theoretical  calcu- 
lations of  Tg  above  model  atmospheres  and  assumptions  of 
linear  increases  in  sea-surface  emissivities  with  wind. 

2.  Actual  analysis  of  Nimbus  6 ESMR  Tg  in  areas  of  known  wind. 

Unfortunately  a calibration  problem  [3]  still  present  in  the  Nimbus  6 
ESMR  (September  1976)  prevented  quantitative  comparisons  of  Tg  and  wind, 
and  the  derivation  of  a sound  empirical  relationship  between  wind  and 
sea-surface  emissivities  needed  to  estimate  winds  from  satellite- 
measured  TB's.  Notwithstanding  the  erroneous  calibration,  analyses  of 
ESMR  Tg  maps  do  show  definite  increases  in  horizontally  and  vertically 
polarized  brightness  temperatures  with  wind.  Such  increases  are  drama- 
tically brought  out  in  a Mistral  occurrence  over  the  Mediterranean  Sea. 

IT]  Sabatini , R.R. , "The  Application  of  the  Nimbus  5 ESMR  to  Sea- 

Surface  Wind  Determination."  EPRF  Technical  Report  6-74  (ESC)  Contract 
No.  N66314-73-C-1572.  Prepared  by  Earth  Satellite  Corporation, 

May.  1974. 

[2]  Sabatini,  R.R. , "Sea-Surface  Wind  Speed  Estimates  from  the  Nimbus  5 
ESMR."  EPRF  Technical  Report  3-75  (ESC)  Contract  No.  N66856-4i20- 
5501.  Prepared  by  Earth  Satellite  Corporation,  February  1975. 

[3]  Communications  with  Dr.  T.  wr‘-''it  the  NASA  ESMR  Experimenter  and 
Dr.  A.  Chang  also  of  NASA  cv  iied  our  suspicions  of  a calibration 
problem.  ESMR  data  are  presently  available  in  calibration  versions 
11,  12,  13,  and  14  all  of  which  are  incorrectly  calibrated.  ESMR 
data  were  first  received  from  NASA  in  versions  11,  12  and  later 

In  version  14  which,  although  gridded  more  accurately  than  the 
previous  versions,  still  contained  a calibration  error. 
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An  error  analysis  on  derived  wind  equations  establishes  the  ac- 
curacy of  wind  speed  determination  from  Nimbus  6 ESMR  under  various 
conditions. 

Section  2 presents  a brief  description  of  the  Nimbus  6 ESMR  experi- 
ment and  data.  In  Section  3 we  discuss  the  surface  and  atmospheric 
effects  on  to  the  microwave  brightness  temperatures  measured  by  the 
Nimbus  6 ESMR,  and  develop  equations  defining  the  relationship  between 
wind  and  Tg.  Section  4 presents  analyses  of  Nimbus  6 ESMR  data.  Section 
5 concludes  the  report  with  a summary  of  results  and  recommendations. 

The  equations  and  model  atmospheres  employed  to  calculate  brightness 
temperatures  above  the  atmosphere  are  presented  in  the  Appendix. 
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2.0  DESCRIPTION  OF  THE  NIMBUS  6 ESMR  EXPERIMENT  AND  DATA 


Only  the  essentials  will  be  presented  here;  for  a more  detailed 
description  of  the  ESMR  experiment  the  reader  should  consult  Section  5 
of  the  Nimbus  6 User's  Guide  [4]. 

The  ESMR  is  one  of  eight  experiments  carried  by  the  Nimbus  6 space- 
craft which  was  launched  in  a near-polar,  circular  (1,100  km),  sun- 
synchronous  (noon  to  midnight)  orbit  on  12  June  1975.  Most  of  the 
experiments,  including  the  ESMR,  are  still  functioning  properly  as  of 
this  writing  (September  1976). 

The  Nimbus  6 ESMR  is  sensitive  to  radiation  in  a 250  MHz  band 
centered  at  37  GHz  (0.81  cm).  The  Nimbus  6 ESMR  measures  both  hori- 
zontally and  vertically  polarized  components  of  the  microwave  radiation 
by  using  two  separate  radiometer  channels.  The  antenna  beam  of  the  ESMR 
scans  ahead  of  the  spacecraft  along  a conical  surface  with  a constant 
angle  of  45°  with  respect  to  the  antenna  axis.  The  beam  scans  in  azi- 
muth ±35°  about  the  forward  direction  in  71  steps,  such  that  beam  posi- 
tion 1 is  35°  azimuth  angle  to  the  right  of  the  spacecraft,  beam  posi- 
tion 36  views  straight  ahead,  and  beam  position  71  views  35°  to  the 
left.  The  instantaneous  f ield-of-view  (half-power  contours)  is  approxi- 
mately an  oval,  20  x 40  km,  nearly  constant  in  size  along  the  scan.  The 
scan  geometry  and  the  tipping  of  the  antenna  5°  forward  of  the  vertical 
axis  of  the  spacecraft  permits  the  antenna  beam  to  intersect  the  earth 
at  a nearly  constant  incidence  angle  throughout  the  scan.  For  the 
expected  spacecraft  pitch  bias  of  +0.6°  the  incidence  angle  varies 
between  an  angle  of  50.8°  for  scan  positions  1 and  71,  and  49.6°  for 
position  36.  The  nearly  constant  incidence  angle  eliminates  the  effects 
of  a varying  nadir  angle  on  sea-surface  emissivity  and  on  atmospheric 
optical  path  lengths,  thus  facilitating  both  qualitative  and  quanti- 
tative data  interpretation. 

The  scanning  geometry  nearly  compensates  for  the  distortions 
introduced  by  the  earth's  curvature.  The  width  of  the  image  area  is 
approximately  1,270  km,  causing  substantial  gaps  in  the  equatorial 
regions  between  successive  orbit  coverages.  These  gaps  decrease  away 
from  the  equator  and  disappear  at  60°  latitude. 

ESMR  data  are  available  to  users  from  the  National  Space  Science 
' Data  Center  (NSSDC)  in  image  format  and  on  Calibrated  Brightness  Temper- 

ature Tapes  (CBTT).  The  ESMR  pictorial  data  and  the  Nimbus  6 catalogs 
[5]  serve  the  purpose  of  determining  coverage  and  selecting  areas  for 
In-depth  analysis  of  brightness  temperature  values  on  the  CBTT. 

All  our  brightness  temperature  analyses  have  been  performed  on 
’ computer  printout  maps  obtained  from  the  NASA-provided  CBT  tapes.  A 

computer  program  in  FORTRAN  and  in  ASSEMBLER  languages  for  the  IBM 
360/168  computer  was  written  to  read  the  CBTT  and  to  map  out  the  cali- 
brated brightness  temperatures  (CBT).  Depending  upon  the  instructions 
inserted  in  a control  card,  CBT  data  to  be  mapped  can  be  selected  on  the 

» [43  Nimbus  Project,  "The  Nimbus  6 User's  Guide,"  NASA,  Goddard  Space 

Flight  Center,  Greenbelt,  Maryland,  1975. 
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basis  of  time  interval  or  latitude-longitude  intervals.  Additionally, 
one  can  choose  to  map  either  the  vertical  or  horizontal  polarization 
brightness  temperatures  (Tn^  or  Tgj^)  or  the  difference  (Tg^  - Tg^). 

The  CBT  maps  contain  all  of  the  original  unaveraqed  values  of  Tg  mapped 
at  a nearly  one-to-one  ratio  between  the  vertical  (along  sub-satellite 
track)  and  the  horizontal  (cross-track)  distances.  The  map  is  at  an 
approximate  scale  of  1:2,200,000.  Gridding  of  the  map  is  achieved  by 
means  of  a rubber-grid  overlay.  The  overlay,  coupled  with  satellite 
subpoint  information  on  the  CBTT  and  printed  out  on  the  CBT  maps, 
greatly  facilitates  gridding  to  nearly  the  accuracy  of  the  resolution  of 
the  ESMR  system.  Landmarks  such  as  coastlines  are  used  to  correctly 
position  the  overlay  grids.  Although  latitude-longitude  information  is 
available  on  the  CBTT  for  every  15th  beam  position  starting  at  position 
six,  these  were  often  found  to  have  errors  of  more  than  one  degree  and 
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3.0  ATMOSPHERIC  AND  SURFACE  EFFECTS  ON  BRIGHTNESS  TEMPERATURES  MEASURED 
?[r37  GHz 

3.1  Introduction 


Theory  and  experiments  have  shown  that  the  microwave  emissivity, 
and  therefore  the  resultant  temperature  of  a sea-surface,  depends 
strongly  on  the  saa  roughness  and  foam  and  spray.  Since  these  are 
related  to  the  wind  velocity  near  the  surface,  sea  brightness 
temperatures  Tg  measured  from  a satellite  can  be  a measure  of  the 
sea-surface  wind  velocity.  The  relationship  between  brightness 
temperature  and  sea-surface  wind  is  complicated  by  factors  which 
can  be  arbitrarily  divided  into  two  groups:  (1)  those  atmospheric 


miVoivyr.  iivjni  OMC 


sea-surface  to  the  satellite  sensor;  and  {l)  those  boundary  layer 
factors  that  influence  the  relationship  between  wind  and  sea  state, 
or  the  effective  emissivity*  of  the  sea-surface. 


Previous  studies  [6,  7]  indicated  that  the  Nimbus  5 ESMR  could 
estimate  sea-surface  winds  to  a five-  to  six-knots  probable  error 
(7.5  to  9 knots  standard  error)  in  clear  areas  and  fully  developed 
open  seas,  given  a cl imatological  estimate  of  atmospheric  water 
vapor  and  sea-surface  temperatures.  These  studies  pointed  out  that 
the  main  difficulty  is  the  definition  of  an  accurate  relationship 
between  wind  and  surface  emissivity,  the  quantity  that  is  actually 
derived  from  Tg  measurements,  and  also  ruled  out  the  use  of  Nimbus 
5 ESMR  Tg  for  wind  estimates  in  rainy  and  cloudy  areas  (except 
cirrus)  unless  the  cloud  liquid  water  could  be  accurately  estimated. 

The  two  polarizations  of  the  Nimbus  6 ESMR  essentially  provide 
two  simultaneous  measurements  of  the  Tg,  which  under  clear  sky 
conditions  permit  estimates  of  the  horizontal  and  vertical  polariza- 
tion sea-surface  emissivities,  which  are  both  related  to  surface 
wind.  Nevertheless,  there  still  remains  the  problem  of  interpre- 
ting the  derived  emissivities  in  terms  of  wind;  even  in  perfectly 
clear  atmospheres  the  ultimate  accuracy  of  the  wind  estimate  relies 
on  how  well  the  sea  wave  spectrum,  foam,  and  spray  (all  of  which 
control  emissivities)  can  be  equated  to  wind.  As  we  shall  show  in 
the  following  sections,  even  with  the  two  Tg  measurements  presence 
of  liquid  clouds  may  cause  unacceptable  errors  in  the  wind  estim- 
ates. 


In  the  following  sections  we  shall  derive  simplified  equations 
for  the  horizontal  and  vertical  brightness  temperatures  at  37  GHz 
In  terms  of  surface  and  atmospheric  parameters,  use  these  equations 
and  an  assumption  of  a linear  relationship  between  wind  and  sea- 
surface  emissivities  to  derive  equations  for  the  wind,  and  ultim- 
ately speculate  on  the  magnitude  of  the  errors  caused  to  the  wind 
estimates  by  errors  in  estimates  of  surface  and  atmospheric  para- 
meters. 

*By  effective  emissivity  is  meant  the  integrated  average  emissivity  of 
an  area  at  least  as  large  as  the  resolution  of  the  sensor,  containing 
an  ensemble  of  waves,  white  caps,  and  spray. 

[6]  Sabatini,  R.R. , 1974,  Op.  Cit. 

[7]  Sabatini,  R.R.,  1975,  Op.  Cit. 
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3 . 2 Derivation  of  Sea-Sur face  Emissivities  and  Wind  Equations 
from  Brightness  Temperatures 


The  brightness  temoerature  above  the  atmosphere  can  be  cal- 
culated for  given  atmospheric  and  surface  conditions  by  means  of 
the  radiative  transfer  equation.  In  our  approach,  called  the  thin- 
atmosphere  approximation  and  detailed  in  the  Appendix,  the  atmos- 
phere is  divided  into  many  layers  each  having  an  average  tempera- 
ture, water  content,  and  absorption  coefficient.  The  radiation 
emerging  from  the  ocean  surface  is  transferred  1 ayer- to- layer  up  to 
the  top  of  the  atmosphere.  In  order  to  derive  a simple  expression 
for  brightness  temperature  in  terms  of  surface  emissivity  and 
atmospheric  and  surface  parameters,  we  solved  the  radiative  trans- 
fer equation  for  clear  and  cloudy  (non-precipitating)  atmospheres. 


aiiu  ciiM:&:>iyiLic^ 
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with  various  sea-surface  temperatures 
fitted  the  results  with  first  order  polynominals . The  resulting 
equations  are  simpler  to  handle  than  the  transfer  equation  and  are 
amenable  to  error  analyses.  Two  sets  of  equations  have  been 
derived:  one  set  for  the  horizontally  polarized  brightness  tem- 
peratures for  which  the  emissivities  were  made  to  vary  from  0.30  to 
0.60,  and  another  for  the  vertical  polarization  brightness  tempera- 
tures for  which  the  emissivities  were  varied  from  0.60  to  0.75. 

The  equations  are  of  the  following  form: 


Tgh  = + kiEhTs  + k2V  + ksC  (3.1) 

^Bv  " '^s^v^s  (3.2) 


where 

Tbv  and  Tg^  are  vertical  and  horizontal  polarization  bright- 
ness temperatures  (°K); 


Ey  and  are  vertical  and  horizontal  polarization  sea-surface 
emissivities  (0.30  1 ^ 0.60),  (0.60  £ Ey  0.75) 

Tj  is  sea-surface  temperature 

V is  atmospheric  water  vapor  in  cm  of  precipitable  water 

C is  non-precipitating  cloud  liquid  water  in  cm  of  precip- 
itable water 

kg  - k3  are  regression  coefficients  for  the  horizontal  polar- 
ization equation 

k4  - k?  are  regression  coefficients  for  the  vertical  polar- 
ization equation 
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Table  3-1  presents  values  of  - k?,  multiple  correlation 
coefficients  Mcc,  and  the  standard  error  of  estimate  SEE  for  a 
variety  of  atmospheres  described  in  the  Appendix.  Figure  3.1  is  a 
plot  of  how  the  brightness  temperatures  Mary  for  the  given  con- 
ditions as  a function  of  emissivity  calculated  with  our  thin  atmos- 
phere assumption,  and  by  means  of  regression  equations  3.1  and  3.2. 
The  plot  gives  an  idea  of  the  goodness  of  fit  of  the  regression 
equations  to  thin  atmosphere  calculations  of  Tg. 

Given  the  brightness  temperatures  and  estimates  of  water  vapor 
and  cloud  liquid  concent,  equations  3.1  and  3.2  can  be  solved  for 
the  sea-surface  emissivities  (E^,,  E^).  The  emissivities  in  turn 
can  be  translated  to  a sea-surface  wind  by  an  appropriate  empirical 
equation  derived  from  Nimbus  6 ESMR.  Unfortunately,  the  inade- 
quately calibrated  ESMR  data  did  not  permit  us  to  derive  reliable 
equations  for  the  sea-surface  emissivities  (both  hurizorital  and 
vertical  polarization)  as  a function  of  observed  winds.  However, 
in  our  subsequent  analyses,  for  the  sake  of  estimating  ranqe  of 
errors  in  the  derived  wind  we  shall  assume  linear  emissivity  versus 
wind  equations  derived  from  37  GHz  observations  from  an  airplane. 

The  vertical  and  horizontal  polarization  sea-surface  emissivi- 
ties are  not  only  determined  by  the  waves,  foam,  and  spray  produced 
by  the  wind,  but  also  by  the  zenith  angle,  and  the  water  tempera- 
tures and  salinity  which  control  the  dielectric  constants.  Figures 
3-2  and  3-3,  derived  from  calculations  presented  by  Paris  [8],  show 
the  sea-surface  emissivities  at  37  GHz  for  calm  conditions  and 
average  salinity  of  32.72  o/oo  as  a function  of  zenith  angle  and 
sea  temperature.  From  these  graphs  we  can  approximate  the  emis- 
sivities for  calm  conditions  at  zenith  = 50“  by: 

Eho  = 1.02059  - 0.162208  x lO’^Tg  - 0.248206  x 10‘^T^  (3.3) 

= 5.4048  - O.O32738T3  + 0.474646  x lO-'^Tj  (3.4) 

where  T^  is  sea-surface  temperature  in  “K.  The  creation  of  waves, 
foam,  and  spray  by  the  wind  increases  these  emissivities.  At  low 
wind  speeds  before  any  foam  forms,  surface  waves  of  dimensions 
comparable  to  the  wavelength  of  observations  are  primarily  respons- 
ible for  the  emissivity  changes.  As  the  wind  speed  increases  so 
does  the  roughness,  but  most  importantly  so  does  the  foam  coverage. 
Graphical  plots  of  observations  of  foam  coverage  with  wind  [9]  show 
a large  scatter,  indicating  that  many  other  factors  affect  foam 
coverage.  These  include  the  duration  and  fetch  of  the  wind,  water 
temperature,  thermal  stability  (difference  in  water  and  air  temper- 
ature), salinity,  and  variations  in  the  surface  tension  of  the 
water  due  to  the  occurrence  of  organic  films.  From  all  available 
observations  of  wind  and  foam  coverage,  Stogryn  [9]  derives  the 
following  equation: 

Tsr  Paris,  J.F. , "Transfer  of  Thermal  Microwave  in  the  Atmosphere."  Vols. 
1,  2.  Dept,  of  Meteorology,  Texas  A&M  University.  Prepared  for  NASA 
NGR-44-001-098,  Office  of  Naval  Research,  Nour  2119  (04),  and 
000  Project  No.  5013,  May,  1971. 

[9]  Stogryn,  A.,  "A  Study  of  Radiometric  Emission  from  a Rough  Sea- 
Surface."  NASA  Contractor  Report  NASA  CR-2038,  July,  1972. 
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TABLE  3.1 


Regression  Coefficients  for  Estimating  Brightness  Temperatures 
at  the  Top  of  the  Atmosphere  with  Equations  3.1  and  3.2 


(1) 

(2) 

(3) 

(4) 

K 

29.4908 

45.7560 

47.4933 

73.9178 

IfO 

^2 

0.8708 

0.7392 

0.7341 

0.5704 

9.8574 

9.2946 

9.2142 

7.0457 

•^3 

0.0 

0.0 

0.0 

624.654 

MCC 

0.9999 

0.996 

0.996 

.992 

SEE 

0.33 

1.80 

2,12 

2.69 

k4 

32.6313 

58.5909 

60.8529 

92.6543 

*^5 

0.8642 

0.7407 

0.7383 

0.5755 

5.2705 

4.7083 

4.7270 

4.0150 

fcc 

0.0 

0.0 

0.0 

351.633 

0.999 

0.995 

0.995 

0.995 

SEE 

0.713 

1.82 

2.01 

1.74 

MCC  - multiple  correlation  coefficient 
SEE  - standard  error  of  estimate  (°K) 

(1)  Sub  Arctic  summer  atmosphere  - clear 
273  < Ts  <283;  0.20  < V < 0.80 

(2)  U.S.  standard  atmosphere  - clear 
280  1 Tj  < 295;  1 .00  ^ V ^ 4.60 

(3)  Mid-latitude  summer  atmosphere  - clear 
290  < T^  300;  1 .50  < V 5.50 

(4)  U.S.  standard  atmosphere  - cloudy 

280  1 Tj  < 295;  1 .00  < V < 4.60;  0.02  < C £ 0.06 


Figure  3-1  Brightness  temperatures  versus  surface  emissivity  at  the 
top  of  a U.S.  Standard  Atmosphere  with  clear  sky,  lOO"', 
clouds  (0.04  cm  precip.  water),  surface  temperature  of 
MO^K,  and  50°  zenith  angle.  Full  lines  are  results  of 
calculations  with  the  thin  atmosphere  approximation  and 
ten  atmospheric  layers  (see  Appendix).  Dashed  lines  represent 
calculations  with  regression  equations  3.1  and  3.2. 
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SEA  SURFACE  EMISSIVITY  AT  37  GH^ 

FOR  CALM  CONDITIONS 
DERIVED  FROM  PARIS  DATA  (1971) 

VERTICAL 

POLARIZ 


Figure  3-3  Horizontal  and  vertical  polarization  calm  sea-surface 
emissivities  at  37  GHz  as  a function  of  temperature, 
derived  from  data  presented  by  Paris  (1971). 
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F = 7.751  X 10-**  W3.231  (3.5) 

where  W is  wind  velocity  in  m/sec,  and  F is  % foam  coverage 
(F  = 100  for  W = 38) 

The  theoretical  determination  of  microwave  emissivity  of  foam 
also  presents  complex  modeling  problems  that  can  best  be  circum- 
vented by  deriving  empirical  equations  from  observations.  Stogryn 

[9]  has  synthesized  all  known  data  on  foam  emissivity  into  a set  of 
empirical  eouations  that  exoress  foam  emissivity  as  a function  of 
temperature,  wavelength,  and  zenith  angle.  At  50'’  zenith  angle  and 
37  GHz  his  equations  reduce  to; 

= 190  (3.6) 

EyTj  =228  (3.7) 

Equations  3.3  and  3.6  show  that  at  = 293°K,  can  vary  from 
0.33  for  calm  conditions  to  0.65  for  100«  foam  cover;  similarly 
from  equations  3.4  and  3.7,  Ey  can  vary  from  0.61  to  0.78. 

Considering  only  the  effects  of  foam  cover  on  the  emissivity, 
equation  3.5  tells  us  that  emissivity  should  increase  slowly  at  low 
to  medium  speed  winds  and  then  much  faster  at  higher  winds.  But  as 
the  wind  increases,  at  low  to  medium  speeds  (^  15  m/sec)  there  is 
also  an  increase  to  the  emissivity  due  to  waves.  The  combination 
of  waves  and  foam  nay  well  be  responsible  for  the  nearly  linear 
effect  of  wind  to  brightness  temperature  observed  with  19.35  GHz 
satellite  microwave  data  by  Sabatini  [10]  in  the  Mediterranean  and 
the  Gulf  of  Tehuantepec  (off  Mexico),  and  observed  with  airborne 
microwave  radiometers  at  various  frequencies  (including  37  GHz)  by 
Webster  et  al.  [11]  in  the  Bering  Sea.  Webster  et  al . [11]  have 
measured  a change  of  0.535  +0.087  °K/meter  per  second  and  1.257  ^ 
0.293  "K/meter  per  second  for  vertical  and  horizontal  polarization 
respectively  at  37  GHz  and  38°  zenith  angles.  Although  these 
changes  were  observed  from  an  aircraft  below  much  of  the  atmosphere 
and  at  a smaller  zenith  angle  than  the  50°  of  the  Nimbus  6 ESMR, 
they  are  indicative  ot  the  magnitude  of  the  change  in  brightness 
temperature  with  wind  expected  in  the  Nimbus  6 ESMR.  Considering 
the  complexity  of  the  emissivity  versus  wind  relationship,  there  is 
no  gain  in  accuracy  in  deviating  from  the  assumption  of  a linear 
Increase  in  brightness  temperature  (as  measured  in  the  absence  of 
an  atmosphere),  and  therefore  a linear  increase  in  emissivity  with 
wind,  at  least  for  winds  up  to  25  m/sec.  With  such  an  assumption 

[10]  Sabatini , R.R.,  1975,  Op.  Cit. 

[11]  Webster,  W.J.  Jr.,  Wilheit,  T.T. , Ross,  D.B.,  and  Per  Gloersen, 
"Analysis  of  the  Convair-990  Passive  Microwave  Observations  of 
the  Sea  States  During  the  Berin  Sea  Experiment,"  Paper  No.  6 in 
"Results  of  the  U.S.  Contribution  to  the  Joint  U.S./USSR  Bering 
Sea  Experiment,"  NASA/GSFC  document  X-910- 74-141 , Greenbelt, 

Maryland,  May  1974. 
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the  vertical  and  horizontal  emissivities  at  37  GHz  can  be  expressed 
by; 


I 


= ^vo  * 

(3.8) 

^h  ■ ^ho  * ^h'^ 

aEv 

AEh 

(3.9) 

W is  wind  in  m/sec;  a^  = aW  and  a^  = aW  are  constants  to  be 
determined  from  observations  of  Tbv  and  in  clear  areas  of  known 
wind,  sea-surface  temperature,  and  atmospheric  water  vapor.  This 
is  achieved  by  plotting  observed  wind  versus  (and  as  cal- 
culated from  equations  3.1  and  3.2.  Unfortunately,  incorrect 
calibration  of  tiimbus  6 ESMR  data  hindered  our  derivation  of  reli- 
able a^  and  ah-  We  can  nevertheless  estimate  a^  and  a^  from  the 
airborne  microwave  observations  of  Webster  et  ai.  [12]  at  37  GHz. 

If  we  assume  that  Tg  = ET^,  their  values  of  0.536  and  1.257  °K/meter 
per  second  for  Alg/AW  roughly  correspond  to  vertical  and  horizontal 
polarization  emissivity  changes  with  wind  (ay  and  a^,)  of  0.002  and 
0.004  per  meter  per  second. 


Equations  for  the  sea-surface  wind  can  now  be  obtained  by 
combining  3.2  and  3.8,  and  3.1  and  3.4; 


^Bh  = *^0  ^ '^i^ho'^’s  '^i^h''’sWh  + k3C 

(3.10) 

Tgy  = ku  + kjEyoT^  + ksa^T^Wy  + k^V  + k7C 

(3.11) 

In  analyzing  equations  3.10  and  3.11  we  founo  that  the  best  pro- 
cedure is  to  obtain  two  independent  estimates  of  the  wind  (W^  and 

Wy)  from  estimates  of  T , V,  and  C.  These  are  then  combined  into 
one  best  estimate  W by  taking  the  weighted  mean  of  W^  and  W^; 

W = [W^/AWh  + W^/AWy]/[l/AW^  + I/AWJ] 

(3.12  a) 

The  total  error  aW  is  then  obtained  from: 

(AW)"  = [(1/AWh)"  + (1/AWy)"]-' 

(3.12b) 

where  and  AWy  are  the  errors  in  and  Wy. 

Let  us  now  briefly  examine  the  accuracy  of  the  above  equations 
for  wind  determination  under  specific  conditions.  We  shall  assume 
that  ay  = 0.002  and  a^  = 0.004,  a U.S.  standard  atmosphere  with 
V = 2 cm  of  precipitable  water,  clear  conditions  or  cloudy  condi- 
tions with  C = 0.04  cm  of  precipitable  water,  and  = 290°K. 

For  clear  atmospheric  conditions  3.12  and  3.13  reduce  to; 

Tgf,  * 118.642  + .857  W,,  + 9.295V  (3.13) 

Tgy  » 191.760  + .430  4.708V  (3.14) 


[12]  Webster  et  al.,  1974,  Op.  Cit. 
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when  we  substitute  the  constants  in  Table  3.1,  and  values  of 
and  Eyo  calculated  from  equations  3.3  and  3.4.  With  errors  of 
1.8°K  in  70^  and  T0v  and  0.1  cm  precipitable  water  in  V,  we  obtain 
errors  of  ± 2.3  m/sec  in  W and  t 4.3  m/sec  in  W^,  yielding  an 
error  of  + 2.0  m/sec  (4  knots)  in  the  average  wind.  If  we  include 
a plausible  * 2°K  error  in  the  sea-surface  temperature,  we  obtain 
an  error  of  * 2.9  m/sec  in  and  i 4.4  m/sec  in  W^,  yielding  an 
error  of  t 2.4  m/sec  (4.8  knots)  in  the  average  wind. 

When  we  consider  the  above  situation  with  cloudy  conditions, 
equations  3.10  and  3.11  reduce  to; 

Tg^  = 130.155  + .662  + 7.046  V + 624.65  C (3.15) 

Tgy  = 196.122  + .334  + 4.015  V + 351.63  C (3.16) 

With  errors  of  2.7°K  in  Tgi^,  1.7°K  in  0.1  cm  precipitable 
water  in  V,  0.01  cm  precipitable  water  in  C,  and  2'’K  in  sea-surface 
temperature,  we  obtain  errors  of  i 10.3  and  -‘11.9  m/sec  in  Wu  and 
Wy,  which  then  combine  to  yield  a total  error  of  7.8  m/sec  (lb 
knots)  in  the  average  wind.  The  error  in  cloud  precipitable  water 
C is  the  biggest  contributor  to  the  total  error;  even  if  we  halve 
this  error  to  0.005  cm  precipitable  water  we  would  still  obtain  a 
total  error  of  4.9  m/sec  (10  knots).  Figure  3.4  suimarizes  the 
errors  in  the  wind  estimates  obtained  from  equations  3.10  and  3.11 
for  clear  and  cloudy  atmospheres  approximating  the  U.S.  standard 
atmosphere. 

The  preceding  error  analysis  assumes  no  errors  in  Tg  measure- 
ments, does  not  include  errors  in  the  emissivity  versus  wind  linear 
relationship  as  expressed  by  equations  3.8  and  3.9,  nor  any  inac- 
curacies introduced  by  the  thin-atmosphere  assumption  and  errors  in 
atmospheric  absorption  coefficients  in  the  derivation  of  3.10 
and  3.11.  All  of  these  factors  would  tend  to  increase  the  wind 
estimate  errors.  Instrument  errors  in  Tg  measurements  can  be 
attenuated  by  averaging  two  or  more  values  of  Tg's;  errors  in  the 
emissivity  versus  wind  relationship  could  be  minimized  by  including 
other  parameters  such  as  fetch  and  air  stability;  finally  the 
errors  of  the  thin-atmosphere  assumption  and  errors  in  the  use  of 
equations  3.1  and  3.2  for  Tg  can  be  eliminated  (especially  in  the 
case  of  a cloudy  atmosphere)  by  solving  the  radiative  transfer 
equation  for  surface  emissivity.  In  the  solution  of  the  radiative 
transfer  equation  we  seek  a surface  emissivity  which  "fits" 
estimated  or  measured  values  of  atmospheric  parameters  and  satellite- 
measured  Tg.  The  determination  of  this  surface  emissivity  is 
therefore  an  iterative  computational  process  which  starts  with  a 
best  estimate  of  E5.  Unless  we  have  very  accurate  measurements  of 
atmospheric  parameters,  such  an  approach  would  yield  no  better 
results  than  the  regression  equation  approach  we  have  used  to 
solve  for  surface  emissivity  (i.e.,  equation  3.1  and  3.2).  With 
our  proposed  approach,  in  a clear  atmosphere  where  the  thin  atmos- 
phere approximation  is  most  accurate  the  wind  estimate  errors  in 
open  seas  (fetch  unlimited  situation)  may  be  kept  to  acceptable 
levels  {*  5 knots)  by  good  estimates  of  atmospheric  water  vapor  (to 
within  S-lO'i  of  total),  sea-surface  temperatures  to  within  2°K,  and 
by  averaging  two  or  more  values  of  Tg's  to  attenuated  instriiment 
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TOTAL  IRNOR  IN  SEA  SURFACE  WIND  (m/Mcl 


Figure  3-4  Errors  in  sea-surface  wind  estimates  using  equations  3.12 
and  3.13  for  a clear  and  cloudy  atmosphere  approximating 
a U.S.  Standard  Atmosphere.  The  lower  horizontal  scale 
represents  the  error  in  water  vapor  content,  AV,  for  a 
clear  atmosphere.  The  upper  horizontal  scale  represents 
the  error  in  cloud  precipitable  water,  aC,  for  a 100% 
cloudy  atmosphere.  An  error  of  - 0.1  cm  of  water  vapor 
content  V is  assumed  for  the  cloudy  atmosphere,  and  of 
± 2®K  in  sea-surface  temperature  for  both  atmospheres. 
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errors.  Such  estimated  wind  accuracies  with  the  Nimbus  6 ESMR  are 
based  only  on  the  analyses  we  just  presented;  verification  must 
await  availability  of  correctly  calibrated  Nimbus  6 ESMR  data. 

In  the  case  of  a cloudy  (liquid)  non-precipitating  atmosphere 
wind  estimate  errors  of  ten  knots  or  greater  may  be  unacceptable. 

In  this  situation  an  iterative  solution  of  the  transfer  equation  to 
obtain  E^  may  be  warranted,  but  only  if  accurate  estimates  of  cloud 
water  content  are  available. 

In  rain  situations  good  sea-surface  wind  estimates  are  prob- 
ably unattainable  because  rain  mostly  obliterates  microwave  radia- 
tion from  the  surface  to  space,  and  therefore  masks  the  effects  of 
surface  emissivity  changes  caused  by  wind.  Figure  3.5  shows  the 
rainfall  absorption  coefficients  at  19.35  GHz  (1.55  cm)  and  37  GHz 
(0.81  cm)  as  calculated  from  equations  given  by  Paris  [13].  At  37 
GHz  rain  absorbs  more  than  twice  than  at  19.35  GHz,  and  even  light 
rain  would  quickly  mask  the  surface  to  the  37  GHz  radiometer. 

The  accuracy  of  sea-surface  wind  determination  in  clear  and 
cloudy  areas  may  be  improved  if  additional  microwave  channel 
measurements  are  available.  Wilheit  et  al.  [14]  have  used  observed 
brightness  temperatures  from  an  airplane  flying  above  the  Bering 
Sea  measured  in  five  microwave  channels  (from  0.81  to  2.3  cm)  to 
derive  estimates  of  the  wind  to  within  1.4  m/sec  even  in  cloudy 
areas.  The  technique  they  employed  involves  the  solution  (in  a 
least  square  sense)  of  five  equations  similar  to  3.10  and  3.11. 

Such  an  approach  from  orbital  heights  will  be- possible  with  a 
satellite  carrying  multichannel  microwave  sensors.  Presently,  to 
test  this  approach  one  could  take  advantage  of  the  Nimbus  5 ESMR  at 
19.35  GHz,  which  when  operative  has  nearly  simultaneous  coverage  of 
the  same  areas  scanned  by  the  Nimbus  6 ESMR.  The  19.35  GHz  Tr 
would  provide  an  additional  equation  of  the  form  of  3.10  and  3.11. 
Assuming  that  such  an  equation  yields  a wind  estimate  with  a 
standard  error  of  about  eight  knots  (four  m/sec),  when  this  wind 
estimate  is  combined  to  the  37  GHz  channel  wind  estimates  with 
standard  errors  of  2.9  and  4.4  m/sec,  we  obtain  a wind  estimate 
with  a standard  error  of  about  two  m/sec  (four  knots).  In  cloudy 
non-precipitating  areas  an  additional  Nimbus  5 ESMR  wind  equation 
having  errors  couipatabie  to  the  two  Nimbus  6 ESMR  equations  (10.3 
and  11.9  m/sec)  would  reduce  the  total  wind  error  to  about  six 
Bi/sec  (12  knots).  The  verification  of  the  accuracy  of  such  a 
three-channel  approach  must  also  await  the  availability  of  cor- 
rectly calibrated  Nimbus  6 ESMR  data. 


tl3]  Paris,  J.F. , 1971,  Op.  Cit. 

[14]  Wilheit,  T.T. , Fowler,  M.G. , Stombach,  G.,  and  Per  Gloersen,  "Microwave 
Radiometric  Determination  of  Atmospheric  Parameters  During  the  Bering 
Sea  Experiment,"  Paper  No.  5 in  "Results  of  the  U.S.  Contribution  to 
the  Joint  U.S. /USSR  Bering  Sea  Experiment,"  NASA/GSFC  document 
X-910-74-141 , Greenbelt,  Maryland,  May  1974. 
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Figure  3-5  Rainfall  absorption  coefficients  as  a function  of  pre- 
cipitation rates  and  rain  content  of  atmosphere  for 
0.81  cm  (37  GHz)  and  1.55  cm  09.35  GHz)  radiation. 
Calculations  were  performed  with  equations  given  by  Paris 
(1971).  Rain  content  and  precipitation  rates  are  related 
by  the  Marshal  1 -Palmer  raindrop  distribution  (1948). 
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4.0  NIMBUS  6 ESMR  OBSERVATIONS  OF  BRIGHTNESS  TEMPERATURES 


A search  for  usable  ESMR  data  over  the  oceans  was  conducted  using 
the  iles  of  ESMR  images  available  at  NASA/GSFC,  Greenbelt,  Maryland, 
and  surface  weather  maos.  We  concentrated  this  first  search  to  the 
North  Atlantic,  an  area  of  dense  wind  observations  and  frequent  cyclon- 
ic storms,  even  in  the  summertime.  The  ESMR  data  selected  from  this 
first  search  covering  only  the  1975  sumer  months  were  finally  received 
from  NASA  at  the  end  of  December  1975,  six  months  after  Nimbus  6 launch. 
Since  these  sunmer  month  data  contained  few  observations  of  high  winds 
above  40  knots,  we  conducted  a second  search  and  put  in  another  order  at 
NASA  for  selected  ESMR  data  over  the  North  Atlantic  during  the  fall  and 
winter  of  1975-1976.  Additionally,  we  ordered  ESMR  data  covering  a few 
possible  Mistral  occurrences  in  the  Mediterranean  sea. 

All  data  were  received  in  the  Calibrated  Brightness  Temperature 
Tape  (CBTT)  format  and  were  mapped  out  with  our  mapping  program.  Maps 
of  Tgy,  Tdu,  and  Tgy  - Tg^  were  produced.  Latitude-longitudes  were 
drawn  on  tne  maps  with  the  aid  of  rubber  grids.  Three  Nimbus  6 daytime 
passes  over  the  North  Atlantic  were  selected  on  the  basis  of  high  winds 
and  good  wind  reports  from  the  first  group  of  ESMR  data. 

Analyses  of  wind  and  brightness  temperatures  in  clear  and  cloudy 
areas  were  conducted  from  the  mapped  Tg  data  and  wind  observations  on 
the  closest  six-hour  surface  maps. 

As  we  proceeded  in  the  analyses  we  noticed  that  observed  brightness 
temperatures  did  not  quite  correspond  to  those  obtained  by  calculations 
with  the  atmospheric  models  described  in  the  Appendix.  The  differences 
were  especially  evident  in  clear  land  and  oceanic  areas  in  which  the 
observed  Tgy  and  Tg^  seemed  to  be  lO-ZO^K  lower  than  the  calculations. 

As  illustrative  examples  of  typical  Tg  values  obtained  from  the  Nimbus  6 
ESMR,  we  shall  present  a few  scans  covering  clear  and  cloudy  areas  as 
shown  by  the  SMS  and  NOAA  satellite  imagery.  Figure  4-1  is  a sector  of 
an  SMS  image  of  the  U.S.  East  Coast  and  Atlantic  areas  on  which  is 
outlined  the  Nimbus  6 ESMR  coverage,  and  two  paths  of  averaged  ESMR 
scans.  Each  swath  represents  the  coverage  of  three  scans,  the  averages 
of  which  are  shown  in  Figures  4-2a  and  4-2b.  The  northern-most  three- 
scan  swath  (Figure  4-2a)  starts  over  land,  crosses  the  Delmarva  Penin- 
sula, skirts  the  New  Jersey  coast,  and  covers  clear  and  cloudy  areas 
over  the  Atlantic.  Tg^  and  Tg^  reach  258°K  and  270°K  respectively  over 
land,  average  124°K  and  192°K  in  clear  ocean  areas  south  of  Long  Island, 
and  reach  up  to  240'’K  and  260°K  over  rain  clouds. 

The  southernmost  swath  (Figure  4-2b),  except  for  some  clouds  in  the 
middle  part,  is  mostly  over  clear  oceanic  areas.  To^  and  Tg  over  clear 
areas  range  from  122-130‘’K  and  196-200°K  respectively  left  or  center, 
and  from  125-145°K  and  196-206°K  right  of  center.  Ocean  brightness 
temperatures  are  higher  than  those  recorded  south  of  Long  Island,  most 
likely  because  of  higher  ocean  temperatures  and  high  atmospheric  water 
content  south  of  the  cold  front.  A gradual  decrease  of  brightness 
temperature  away  from  the  center  of  the  swath  is  apparent,  especially 
In  the  Tgh- 
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Figure  4-2  Average  horizontal  and  vertical  brightness  temperatures 
of  three  scan  lines  occurring  over  the  U.S.  East  Coast 
and  Atlantic  area  at  about  40°N  (a),  and  at  about  32°N 
(b).  The  position  of  the  two  sets  of  scan  lines  is  also 
shown  on  the  SMS-1  image  in  Figure  4-1. 
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For  comparison  to  the  observed  Tp's  we  call  attention  to  Figure  3-1 
which  shows  the  results  of  Tp  calculations  for  a U.S.  standard  atmos- 
phere and  sea-surface  temperature  of  290°K  obtained  by  our  thin  atmos- 
phere approximation  (see  Appendix).  These  calculations  yield  clear 
areas  Tg^  and  values  of  about  204°K  and  144'’K  at  vertical  and  hori- 
zontal polarization  sea-surface  emissivities  of  0.61  and  0.34.  These 
are  12°  and  20°K  higher  than  the  comparable  Oimbus  6 ESMR  averages  of 
192°K  and  124°K  south  of  Long  Island.  Other  data  checks  with  calcu- 
lat  ions  made  with  our  atmospheric  model,  comparisons  of  our  model 
results  to  calculations  performed  by  Stogryn  [15],  and  to  calculations 
made  with  formulas  given  by  Paris  [16]  at  55°  Nadir,  also  indicated  that 
the  Nimbus  6 ESMR  produces  Tg^  and  Tg^  10-20°K  lower  than  what  they 
should  be.  Subsequent  conversations  with  NASA  officials  at  GSFC,  who  at 
about  the  same  time  (May  1976)  were  experiencing  the  same  difficulty 
with  the  data,  confirmed  our  suspicions.  According  to  Dr.  Wilheit  [17] 
the  ESMR  experimenter,  the  problem  is  due  to  the  heating  of  the  sun 
shining  on  the  microwave  antenna  causing  an  erroneous  calibration  of  the 
instrument.  NASA  is  presently  working  on  deriving  a corrected  cali- 
bration which  would  eliminate  this  heating  effect.  An  analysis  we 
conducted  with  some  of  the  data  shows  that  this  calibration  error  is 
more  pronounced  in  the  horizontal  polarization,  and  is  a function  of 
beam  position  (which  determines  antenna  aspect  position  relative  to 
sun).  Figure  4-3  presents  Tgy  and  Tgv^  values  averaged  for  each  beam 
position  for  swaths  spanning  the  south  and  north  Atlantic  Ocean.  To 
eliminate  thick  clouds  and  rain  areas  which  would  introduce  gross  devi- 
ations in  these  limited  area  averages,  we  included  only  Tg^  below  170°< 
and  Tgy  below  210"K.  The  four  swaths,  two  in  the  daytime  (24  November 
1975,  29  February  1976)  and  two  in  the  nighttime  (14  December  1975, 

21  February  1976)  were  chosen  at  random  from  the  ESMR  data  we  had 
ordered  from  NASA.  All  averages  show  a peak  near  the  center  of  the 
swath,  a steep  decrease  away  from  the  center,  and  a small  rise  near  the 
edges.  For  this  limited  data  sample,  the  peak-to-trough  variations  are 
about  7-ll°K  for  Tg^  and  12-17°K  for  Tg^,  with  the  data  near  the  center 
of  the  scan  containing  the  least  relative  errors.  The  corrections  to  be 
applied  are  obviously  a function  of  antenna  position  relative  to  the 
sun,  which  in  turn  is  determined  by  beam  position,  latitude,  and  time  of 
year.  No  attempt  was  made  in  our  work  to  correct  for  this  calibration 
error. 

Calibration  errors  render  impossible  the  use  of  the  Nimbus  6 ESMR 
data  for  deriving  valid  relationships  between  sea-surface  winds  and 
brightness  temperatures.  Nevertheless,  as  a preliminary  effort  to 
evaluate  the  data  at  hand  for  sea-surface  wind  estimation  we  have  con- 
sidered the  analysis  of  a Mistral  occurrence  in  the  Mediterranean  sea 


[16] 

[17] 


Stogryn,  A.,  1972  Op.  Cit. 

Paris,  J.F.,  1971,  Op.  Cit. 
Communication  with  Dr.  T.  Wilheit,  May, 


1976. 
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Figure  4-3  Average  brightness  temperatures  as  a function  of  beam 

position  for  four  Nimbus  6 ESMR  swaths  spanning  the  south 
and  north  Atlantic  ocean.  The  four  swaths,  two  in  the 
daytime  (24  November  1975,  29  February  1976)  and  two  in 
the  nighttime  (14  December  1975,  21  February  1976)  were 
chosen  at  random  from  ESfIR  data  furnished  by  NASA.  The 
solid  and  dashed  lines  represent  1-2-1  smoothing  of  the 
actual  averages. 
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and  the  analysis  of  large  area  averages  of  wind  versus  Tg.  The  large 
area  analysis  would  attenuate  the  calibration  errors. 

We  shall  first  present  a Nimbus  6 ESMR  coverage  of  a Mistral  occur- 
rence on  15  January,  1976,  because  it  is  probably  the  most  convincing 
evidence  we  have  that  Nimbus  6 ESMR  can  be  used  for  sea-surface  wind 
determination.  Figure  4-4  and  4-5  show  the  gridded  and  maps 
obtained  from  a Nimbus  6 pass  over  the  Mediterranean  on  15  January 
1976,  between  10:31  and  10:33  GMT.  Recorded  on  the  maps  are  the  reported 
winds  at  1200  GMT  which  range  up  to  45  knots.  Figure  4-6  is  a NOAA 
satellite  image  of  the  area  taken  at  8:48  GMT,  on  which  is  outlined  the 
ESMR  coverage.  The  surface  weather  map  for  12  GMT,  15  January  (Figure 
4-7)  contains  the  typical  Mistral  synoptic  conditions  with  a large 
anticyclone  centered  over  western  Europe  and  the  Atlantic  and  a trough 
crossing  the  Italian  peninsula.  Both  the  T^^  and  the  T maps  show  an 
area  of  relatively  higher  temperatures  in  tKe  Gulf  of  LiJ^n,  yet  the  NOAA 
satellite  image  and  the  1200  GMT  surface  map  show  little  or  no  clouds  in 
the  area.  Available  ship  reports  at  1200  GMT  show  one-  to  two-tenths 
cirrus  and  one-tenth  cumulus,  which  would  hardly  affect  the  brightness 
temperatures. 

Relatively  higher  brightness  temperatures  in  the  Gulf  of  Lion  area 
during  Mistral  occurrences  have  also  been  reported  in  a previous  study 
with  Nimbus  5 ESMR  [18],  and  are  thought  to  be  caused  by  the  higher 
winds  of  the  Mistral  funnel ing  into  the  Gulf  of  Lion  and  causing  rougher 
seas  in  this  area.  Ships  in  and  around  the  area  of  the  Mistral  reported 
winds  of  35  to  45  knots.  In  one  instance,  near  the  coast  of  France  at 
about  41. 5N  and  3.4E,  the  high  wind  report  is  not  reflected  by  higher 
Tb’s.  One  may  think  that  calibration  errors  and  the  88  minutes  differ- 
ence between  the  wind  observation  and  the  Tg  map  could  be  the  cause  of 
such  discrepancy.  A more  likely  explanation  is  the  nearness  of  land,  a 
reminder  that  what  the  ESMR  "sees"  is  only  the  effect  of  the  wind  upon 
the  waters,  i.e.  waves,  foam,  and  spray.  Near  the  coastline  where  the 
fetch  is  short,  the  effect  of  the  Mistral  wind  upon  the  waters  is  not  as 
pronounced  as  in  the  open  seas  and  therefore  Tg's  are  not  as  high. 

Figure  4-8  represents  the  average  Tg^  and  Tg^  of  three  scans  across 
the  Mistral  area  (see  Figure  4-7).  These  averages  also  bring  out  the 
effect  of  the  Mistral  in  the  Gulf  of  Lion,  which  is  evidenced  by  nearly 
25°K  and  20°K  increases  in  the  Tg^  and  Tg^  averages  respectively. 
Partially  masking  and  distorting  the  Mistral  effect  is  of  course  the 
calibration  error.  If  we  can  assume  that  the  curves  shown  in  Figure  4-3 
are  representative  of  the  relative  calibration  error  across  beam  posi- 
tions 1 to  71 , then  the  calibration  error  would  emphasize  the  Mistral 
effect  near  the  left  edge  (beams  65-53)  by  steepening  the  west-east  Tg 
gradient  here  by  about  4-6°K;  and  would  de-emphasize  the  effect  west  of 
the  apparent  Tg  maximum  (beam  53)  by  lessening  the  east-west  Tg  gradient 
by  about  4-5°^  in  the  following  12  beam  positions  (52-41). 

TTeT  Sabatini,  R.R.,  1974,  Op.  Cit. 
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Figure  4-7  Surface  weather  map  for  1200  GMT  15  January  1976. 
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Figure  4-8  Average  horizontal  and  vertical  briohtness  temperatures 

of  three  scan  lines  crossing  the  Mistral  area  on  15  January 
1976.  The  position  of  the  scan  lines  is  also  shown  on 
the  NOAA  satellite  image  in  Figure  4-6. 
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Because  of  calibration  errors  and  insufficient  wind  data,  no  quan- 
titative evaluation  of  Tg  versus  wind  is  possible  from  this  single 
Mistral  event.  Nevertheless,  this  Mistral  occurrence  does  show  a pro- 
nounced effect  on  37  GHz  brightness  temperatures  similar  to  the  effect 
observed  on  the  19.35  GHz  Tg's  in  other  Mistral  occurrences,  and  indi- 
cates that  further  analyses  with  correctly  calibrated  37  GHz  Tg  data  may 
yield  useful  statistical  relationships  between  wind  and  Tg's. 

The  next  analysis  involves  large  area  comparisons  of  wind  and 
brightness  temperatures. 

Three  Nimbus  6 TSMR  daytime  passes  over  the  North  Atlantic  were 
selected  for  this  analysis  from  the  first  batch  of  ESMR  data  received 
from  NASA  at  the  end  of  December  1975.  The  selection  was  made  on  the 
basis  of  high  winds  and  good  wind  reports.  The  passes  occurred  on  5 
August,  11  September,  and  13  September  1975.  Maps  of  Tn^,  Tr;^,  and 
Tgy  - Tg^  covering  the  North  Atlantic  area  were  produced  wit^  our  map- 
ping program  and  gridded.  The  ESMR  data  which  v/ere  mapped  occurred  from 
about  12:20  to  12:30  GMT  on  5 August,  from  about  13:A5  to  13:55  GMT  on 
11  September,  and  from  about  15:45  to  15:55  GMT  on  13  September. 

Surface  winds  over  the  ocean  were  estimated  from  observations  and  from 
pressure  gradients  by  means  of  the  geostrophic  approximation  corrected 
for  surface  friction.  The  12:00  GMT  surface  maps  and  wind  observations 
were  used  for  the  August  5 and  September  11  cases.  Both  the  12:00  GMT 
and  the  18:00  GMT  maps  v/ere  used  to  interpolate  the  position  of  the 
Isobars  for  the  time  of  the  September  13  ESMR  data  (about  16:00  GMT). 

Isotachs  were  drawn  at  five-  and  ten-knot  intervals  and  transferred 
to  the  gridded  Tg  maps.  Average  Tg^  - Tn^  were  then  calculated  for  the 
area  enclosed  by  each  ten-knot  interval  (or  five-knot  interval  where 
wind  data  warranted  it).  F'gure  4-9  shows  a plot  of  the  results.  Each 
point  is  an  average  Tg^  - Tg^  plotted  against  the  middle  of  the  wind 
interval;  i.e.,  the  average  of  all  Tg^  - Tg^  falling  within  0-to  10- 
knot  isota  's  v/ere  plotted  at  5 knots,  etc.  The  graph  in  Figure  4-9 
shov/s  the  ..pected  decrease  in  Tg^  - Tg^  as  the  wind  increases.  How- 
ever, great  caution  should  be  used  in  interpreting  these  results.  First 
of  all,  no  allowances  were  made  for  the  calibration  error  discovered 
after  this  analysis  was  completed;  second  and  most  importantly,  the 
averages  include  data  over  all  clear,  cloudy,  and  rain  conditions.  The 
correspondence  between  Tg,/  - Tg^  and  wind,  if  any,  may  therefore  be 
hidden  by  the  fact  that  there  is  an  increase  in  cloudiness  and  pre- 
cipitation in  areas  where  one  also  expects  an  increase  in  wind  (toward 
the  center  of  a low  pressure  system).  Clouds  and  rain  lessen  the  con- 
tribution of  the  surface  polarized  radiation,  and  therefore  tend  to 
attenuate  Tg^  - Tg^,. 

Figure  4-10  a.b  shows  the  average  Tgv  and  Tg^  plotted  against  the 
middle  of  the  wind  interval  for  the  August  5 case  only.  Since  the 
calibration  error  was  discovered  at  about  the  time  of  this  analysis,  we 
decided  not  to  proceed  with  further  analyses  on  the  other  two  cases. 
Again,  although  both  Tgy  and  Tgg  show  the  correct  trend  with  wind,  we 
cannot  come  to  any  firm  conclusion  because  of  the  inclusion  of  cloud 
areas.  A valid  analysis  of  Tg  versus  wind  must  await  correct  cali- 
bration of  the  data,  and  should  exclude  cloud-covered  areas.  For  such 
an  analysis,  clear  areas  can  best  be  determined  by  temperatures  obtained 
from  the  Temperature  Humidity  Infrared  Radiometer  (THIR),  also  flying  on 

the  Nimbus  6 satellite. 
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Figure  4-9  Plot  of  large  area  averages  of  Tg^  - Tg^  as  a function  of 
wind.  The  ESMR  data  used  in  this  analysis  occurred  over 
the  North  Atlantic  from  about  12:20  to  12:30  GMT  on 
5 August  1975,  from  about  13:45  to  13:55  GMT  on  11  Sep- 
tember 1975,  and  from  about  15:45  to  15:55  GMT  on  13  Sep- 
tember 1975. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  presence  of  a calibration  error  has  prevented  a quantitative 
evaluation  of  the  applicability  of  the  Nimbus  6 ESMR  to  sea-surface  wind 
determination.  Comparisons  of  Ip's  produced  by  the  Nimbus  6 ESMR  to 
Tg's  calculated  over  model  atmospheres  show  that  calibration  errors 
cause  Tgy  and  ^Bh  to  be  10-20  °K  lower  than  what  they  should  be.  The 
error  is  determined  by  the  sun-antenna  anqle,  and  therefore  shows  up  as 
a function  of  antenna  beam  position  with  the  error  increasing  toward  the 
horizons. 

A preliminary  error  analysis  on  wind  equations  which  assume  a 
linear  relationship  between  wind  and  sea-surface  emissivities  shows  that 
in  open,  wel 1 -developed  seas,  surface  winds  less  than  50  knots  may  be 
estimated  from  Nimbus  6 ESMR  to  an  accuracy  of  about  5 knots  in  clear 
areas  and  10  to  16  knots  in  cloudy  non-precipitating  areas.  The  addi- 
tion of  Nimbus  6 ESMR  data  would  improve  the  wind  estimates  to  an 
accuracy  of  4 knots  in  clear  areas  and  7 to  1 1 knots  in  cloudy  non- 
precipitatinq  areas.  These  accuracies  can  further  be  improved  by 
solving  for  surface  emissivities  by  iterative  solutions  of  the  microwave 
radiative  transfer  equation;  this  method  is  warranted  only  if  atmos- 
pheric parameters  are  accurately  known. 

The  high  sensitivity  of  the  Nimbus  6 ESMR  to  rain  precludes  its  use 
for  sea-surface  wind  estimation  in  rainy  areas. 

The  analysis  of  Nimbus  6 ESMR  Tg’s  during  a Mistral  occurrence 
shows  that  even  in  the  presence  of  a calibration  error  the  Nimbus  6 ESMR 
can  qualitatively  detect  the  effect  of  the  wind  upon  the  sea.  This 
effect  causes  relatively  hiqh  horizontal  and  vertical  brightness  temp- 
eratures in  areas  of  expected  higher  winds.  Large  area  analysis  of  Tg's 
and  wind  also  show  increases  in  Tgy  and  Tg^  with  wind,  and  a decrease  in 
Tgv  ■ Tg^.  This  large  area  analysis  is  nevertheless  inconclusive, 
because  of  the  natural  presence  of  an  increase  in  cloud  cover  and  rain 
with  winds,  which  also  increase  Tg's. 

It  is  concluded  that  Tg's  versus  winds  analyses  should  be  pursued 
further  when  correctly  calibrated  Nimbus  6 ESMR  data  are  available. 
Concurrent  Nimbus  5 ESMR  data  should  also  be  used  to  derive  and  test  an 
improved  set  of  three  wind  estimation  equations. 

Analyses  for  deriving  sea-surface  emissivity  versus  wind  equations 
should  concentrate  on  clear  areas  in  which  atmospheric  effects  can  be 
best  corrected  for.  Advantage  should  be  taken  of  simultaneous  Nimbus 
THIR  measurements  for  determining  clear  areas. 

The  wind  equations  derived  by  the  proposed  method  involving  best- 
fit  equations  of  Tg  versus  surface  and  atir.ospheric  parameters  should  be 
tested  against  a second  method  employing  iterative  solutions  of  the 
microwave  radiative  transfer  equation  for  surface  emissivity.  Tests 
should  be  made  over  clear  and  cloudy  non-precipitating  areas  in  which 
surface  winds  and  radiosonde  ascents  are  available. 
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APPENDIX 


EQUATIONS  FOR  BRIGHTNESS  TEMPERATURE  CALCULATIONS 
AND  DlSCR'IPTT'QN  O'F'  MODEL  ATMOSPHERES  ^ ' 


^ ^ The  T_r ansfer  Equation  for  Microwave  - Th e_ T hj_n_ A tmqsj)ji e rq^i - 
mation 

The  microwave  brightness  temperature  Tg  measured  by  a radiometer  on 
a satellite  may  be  expressed  as; 

Tg  = T(A)  + T(S)  + T(R)  (1) 

where  T(A)  represents  the  upward  emission  from  the  atmosphere,  T(S)  is 
the  attenuated  surface  emission,  T(R)  is  the  attenuated  portion  of  the 
reflected  downv/ard  emission,  T{D).  These  four  terms  may  be  expressed 
as: 


T(A)  = I"  T(2)  a(2)exp  ( - a(z)  d2)  d2 

T(S)  = CgTs  exp  ( - /“  n.(2)dz) 

T(R)  = (1  - Es)  exp  ( - at(2)d2)  T(D) 

and  T(D)  = /"  (T)2  a (2)  exp  ( - /”  a(z)dz)  dz 

In  the  above  equations  and  T^  are  surface  emissivity  and  temperature, 
T(2)  and  1(2)  are  atmospheric  temperature  and  absorption  coefficient  as 
a function  of  height,  2. 

The  terms  of  equation  (1)  can  be  evaluated  numerically  by  dividing 
the  atmosphere  into  a finite  number  of  layers  of  thickness  tz,  each  with 
an  average  temperature  and  absorption  coefficient.  For  our  calculation 
we  used  a ten-layer  model. 

A . 2 Absorption  Coefficients 

Water  in  all  its  forms  and  molecular  oxygen  are  the  main  absorbers 
of  microwave  radiation  in  the  atmosphere.  The  microwave  absorption 
coefficient  of  a layer  is  the  sum  total  of  the  absorption  coefficients 
of  the  constituents: 

a = + ai  + o.  + 

where  oq  is  the  absorption  coefficient  of  molecular  oxygen  and  a , , 

a-j,  are  the  absorption  coefficients  of  water  in  vapor,  small  Yiquid 
drops,  ice,  and  rain  forms;  ci^,  is  the  absorption  coefficient  of  small 
melting  ice  spheres.  By  "small"  are  intended  cloud  particles  with  an 
average  radius  smaller  than  50  microns  for  which  the  Rayleigh  absorption 
coefficient  does  not  apply.  The  Rayleigh  approximation  does  not  apply 
to  the  larger  particles  constituting  rainfall.  The  absorption  coeffici- 
ent of  rain  has  to  be  calculated  by  means  of  the  complex  Mie  theory. 
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A . 2 . 1 Absorption  Coe fficient  of  Water  Vapor 


The  water  vapor  absorption  is  calculated  by  equations  given  by 
Westwater  [1,  2].  These  equations  are: 

ajp  = (318/T)2.5  exp  (-644/T  + 644/318)o?Ci 

{L/[{v  - 0)2  + l2)]  + L/[(0  + 0^)2  + l2)]} 

+ (318/T)C°o2l 

where  = water  vapor  absorption 

p = water  vapor  density  (gm'^) 

L = (Av/C^)w  = line  width  (cm"*) 

(o,v^)  = wave  number  and  resonant  wave  number  (cm~^) 

T = absolute  temperature  (K) 

The  pressure,  P(mb)  and  temperature  dependences  of  the  line  width  are 
expressed  by: 

L = (P/1013. 25)(318/T)-625  a(i  + bp) 

Numerical  values  of  the  constants  are; 

C 

Cl  = .0008312 
C2  = .01402 
a = .08478 
b = .00708 
Oq  = .7417  cm"^ 

A . 2 . 2 Absorption  Coefficient  of Non-Precipitating  Clouds 

The  absorption  (or  extinction)  coefficient  for  drop  size  distri- 
butions typical  of  non-precipitating  clouds  can  be  expressed  in  terms  of 
the  liquid  water  content  of  the  cloud  particles  as  shown  by  Deirmendjian 
[3].  For  wavelengths  much  larger  than  the  average  drop  size  in  a cloud, 
scattering  is  negligible  and  the  attenuation  can  be  approximated  by 
Rayleigh's  theory.  The  Rayleigh  model  of  absorption  is  applicable  to 
clouds  whose  droplets  range  from  a few  microns  to  a few  hundred  microns 
and  the  absorption  coefficient  is  given  by; 

a = 1.885Q/A  H-K)  (2) 

where  a = absorption  coefficient  (km"i) 

X = wavelength  (cm) 

Q = liquid  water  content  (g/m^) 

and  K = (e  - l)/(c  + 2)  (3) 

where  e = complex  dielectric  constant.  The  notation  If-K)  means  the 
negative  of  the  imaginary  part  of  the  complex  variable  K. 

The  dielectric  constant  of  water  is  given  as  a function  of  wave- 
length and  temperature  T(°K),  by  Grant,  et  al.  [4]  as: 
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(4) 


e * + (e^  - e„)/[l  + (iX^/X)®] 

where  c„  = 4.5 

= 32155. 45/T  - 29.62 
logx  = 921.0936/T  - 2.9014 
I = 0.98 
1 = 

Solving  for  the  imaginary  part  of  K by  introducing  (4)  into  equation 
(3)  and  using  the  identity; 

i = exp  (i;T/2)  = cos  ti/2  + i sin  ii/2 

we  obtain  for  water  particles: 

I{-K}  = 3£^  A^sin  (ttS/2)/D 

and  A = x,/x 

D = + {e„  + 2)2n  + a2^) 

+ 2(g„  + 2)2  A?'cos(rp/2) 

+ 2g^(g„  + 2)  [1  + ABcos(ne/2)] 

The  dielectric  constant  of  ice  is  given  by  Deirmendjian  [3]  as: 

G = 3.1684  - .0085441 
and  therefore, 

H-K}  = .0009595 

and, 

= .0023  lan‘1  at  37  GHz 

For  small  drops  made  of  ice  and  water,  K is  given  by  Van  de  Hulst 
[5],  as; 

(g2  - 1)(gi  + 2e2)  + qM2e2  + l)(ei  - £2) 

K(G],C2tq)  “ 

(cj  + 2)(gi  + 2g2)  + q^(2G2  - 2)(gi  - E2) 

%<here 

Cj  * inner  dielectric  constant 
E2  = outer  dielectric  constant 
q = ratio  of  inner  to  total  radius 
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A . 2 . 3 Absorption  Coefficient  of  Rain 

When  the  particle  size  approaches  the  wavelength  of  the  incident 
radiation  scattering  is  appreciable  and  must  be  considered  along  with 
the  absorption,  and  therefore  the  attenuation  cannot  be  approximated  by 
the  Rayleigh  theory.  The  more  complicated  Mie  model  of  absorption  and 
scattering  must  be  used.  Paris  [6]  used  the  Mie  theory  to  calculate  the 
volume  absorption  coefficient  of  a pol /dispersive  cloud  of  hydrometeors 
obeying  the  Marshal  1 -Palmer  drop-size  distribution  [7].  The  Marshall- 
Palmer  drop-size  distribution  is  fairly  representative  of  the  drop-size 
distribution  found  in  most  sub-tropical  clouds  and  can  be  readily  ex- 
pressed in  a rate  of  precipitation  (R,  mm/hr)  or  a liquid  water  content 
(M,  g/m^)  by  the  relationship  M = Paris  fitted  his  results 

to  a regression  equation  of  the  form; 

where  the  absorption  coefficient  of  rain  n (m"M  is  for  a Marshall- 
Palmer  raindrop  distribution,  M is  the  content  of  liquid  water  (gm"^)  and 
k and  d are  constants  which  are  a function  of  temperature,  T(C).  At 
37  GHz  these  are  [6]: 

k = 0.74322  X lO'^  - 0.33730  x lO'^T  + 0.60416  x lO'^T^ 

- 0.35503  X 10-^3  + 0.63349  x lO'^n** 

d = 1.0318  - 0.54367  x lO'^T  - 0.19236  x 10“^^ 

- 0.66322  X 10-^T3  + 0.22333  x lO'^T** 

A . 2 . 4 Absorption  Coefficient  of  Molecular  Oxygen 

The  absorption  coefficient  for  molecular  oxygen  a (m"^)  is 
calculated  by  equations  given  by  Meeks  and  Lilley  [8].  tn  meteorologi- 
cal units: 

a = 4.6182  X 10'13  (p  v^/T^) 
o 

• t exp[-2. 06844  N (N  + 1)/T], 

N = 1 
odd 
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where 


S 

N 


r 2(N2  + N + L)  (2N  + 1)  c N(2N  + 3) 
"o  ^ ■ 

N(N+1)  N + 1 


-f  Fn- 

N 


p is  the  atmospheric  pressure  in  mb,  T is  atmospheric  temperature  in  °K, 
and  V is  frequency  of  radiation  in  Hz. 

F^,  F are  pressure  broadening  effects  which  are  of  the  following 
forms;  ’ 


and 


Ni  ■ 


Hi  - v)^  + 


Ao 


“0 

+ V ) '■  + 


(N  = 1,  3,  5,  45)  are  the  resonant  frequencies  (Hz)  for 

molecular  oxygen  listed  in  Table  A-1,  u is  the  frequency  (Hz)  of  the 
radiation  under  consideration,  and  Aq  is  the  line-width  parameter  (Hz) 
for  molecular  oxygen  given  by  Meeks  and  Li  1 ley  as  follows; 


Aq  = 1.4625  X 10^  p (300/T)0.85  (o.2i  + o. 78  f) 


where  p is  the  atmospheric  pressure  (mb),  T is  the  atmospheric  temper- 
ature (K),  and  f is  a factor  that  expresses  the  relative  effectiveness 
of  nitrogen-oxygen  collisions  as  compared  to  oxygen-oxygen  collisions  in 
producing  broadening.  The  value  of  f changes  from  0.75  for  very  low 
pressure  to  0.25  for  high  pressures.  Meeks  and  Lilley  (1963)  give  the 
following  empirically-derived  form  for  f; 


f = 


^0.25;  p > 356  mb 
■ 0.25  + 0.435  (2.551 
^0.75:  p 25.  3 mb 


logjpp);  25. 3 < p < 356  mb 


We  now  have  the  complete  set  of  equations  required  to  calculate  the 
brightness  temperature  through  the  atmosphere  by  the  thi n-atmosphere 
method. 


A. 3 Atmospheric  Models 

The  atmospheric  models  used  for  our  brightness  temperature  cal- 
culations are  divided  into  10  layers  from  the  surface  to  15  Km.  The 
layers  are,  starting  from  the  surface;  1 , 1 , 1 , 1 , 1 , 1 ,2,2,2,3  Km  thick, 
each  characteri zed  by  an  average  temperature,  pressure,  and  water  vapor. 
Three  types  of  atmospheres  were  used  as  the  basis  for  deriving  regres- 
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Sion  equations  for  Tp.  These  are  (1)  the  U.S.  standard  atmosphere 
(1962),  (2)  a sub-arctic  summer  atmosphere,  and  (3)  a midlatitude  summer 
atmosphere.  All  three  model  atmosphere  were  taken  from  Valley  [9]. 
Clouds  were  inserted  only  in  the  U.S.  standard  atmosphere  from  one  Km  to 
near  the  freezing  level . 

To  obtain  more  atmospheric  models,  the  vapor  and  cloud  contents 
were  varied  from  approximately  50S  to  200'  higher  than  the  values  shown 
in  Tables  A-2  to  A-4.  Atmospheric  temperatures  remained  constant. 
Calculations  were  performed  for  sea-surface  emissivities  varying  from 
0.30  to  0.80  and  sea-surface  temperatures  varying  from  273°K  to  300‘K. 


A-6 


TABLE  A-1 


Resonant  frequencies  for  the  absorption  of  microwaves  by  molecular 
oxygen  (after  Meeks  and  Lilley,  1963) 


''n+ 

N 

''N+ 

'"N- 

N 

(GHz) 

(GHz) 

(GHz) 

(GHz) 

1 

56.2648 

118. 7505 

25 

65.7626 

53.5960 

3 

58.4466 

62.4863 

27 

66.2978 

53.0695 

5 

59.5910 

60.3061 

29 

66.8313 

52.5458 

7 

60.4348 

59.1642 

31 

67.3627 

52.0259 

9 

61 .1506 

58.3239 

33 

67.8923 

51 .5091 

11 

61 .8002 

57.6125 

35 

68.4205 

50.9949 

13 

62.4112 

56.9682 

37 

68.9478 

50.4830 

15 

62.9980 

56.3634 

39 

69.4741 

49.9730 

17 

63.5685 

55.7839 

41 

70.0000 

49.4648 

19 

64.1272 

55.2214 

43 

70.5249 

48.9582 

21 

64.8779 

54.6728 

45 

71.0497 

48.4530 

23 

65.2240 

54.1294 
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TABLE  A-2 


U.S. 

STANDARD  ATMOSPHERE 

FOR  Tg 

CALCULATIONS 

Layer 

Approx . 
Height 
(km) 

Pressure 

(mb) 

Temp. 

m 

Water 

Vapor 

(gm-3) 

Cloud 

Liquid  Water 
(gm-3) 

1 

0.5 

956 

285 

10.56 

2 

1.5 

847 

278 

6.73 

0.2 

3 

2.5 

748 

272 

4.48 

0.2 

4 

3.5 

659 

265 

2.15 

5 

4.5 

579 

259 

1.33 

6 

5.5 

506 

252 

0.95 

7 

7.0 

411 

243 

0.37 

8 

9.0 

308 

230 

.80 

9 

11.0 

227 

217 

.01 

10 

13.5 

154 

217 

.00 

Total 

Precipitable  Water  (cm) 

2.71 

0.04 

A-8 


TABLE  A-3 


MIDLATITUDE  SUMMER  ATMOSPHERE  FOR  Tg  CALCULATIONS 


Layer 

Approx. 

Height 

(Km) 

Pressure 

(mb) 

Temp. 

(°K) 

Water 

Vapor 

(gm-3) 

1 

0.5 

958 

289 

13.51 

2 

1.5 

852 

285 

10.56 

3 

2.5 

756 

280 

7.68 

4 

3.5 

669 

275 

5.50 

5 

4.5 

591 

269 

3.62 

6 

5.5 

521 

261 

0.80 

7 

7.0 

426 

248 

0.40 

8 

9.0 

324 

232 

0.10 

9 

11.0 

243 

217 

0.02 

10 

13.5 

166 

217 

0.00 

Total  Precipitable  Water  (cm) 

4.27 
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TABLE  A-4 


SUBARCTIC  SUMMER  ATMOSPHERE  FOR  Tg  CALCULATIONS 


Layer 

Approx. 

Height 

(km) 

Pressure 

(mb) 

Temp. 

{°K) 

Water 

Vapor 

(gm-3) 

1 

0.5 

958 

285 

7.50 

2 

1.5 

845 

279 

5.10 

3 

2.5 

747 

274 

3.45 

4 

3.5 

658 

269 

2.20 

5 

4.5 

579 

263 

1.35 

6 

5.5 

507 

257 

0.77 

7 

7.0 

413 

246 

0.28 

8 

9.0 

311 

232 

0.02 

9 

11.0 

230 

225 

0.01 

10 

13.5 

158 

225 

0.00 

Total  Precipitable  Water  (cm) 

2.10 
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